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ABSTRACT 

We present spectrally resolved observations of the young multiple system T Tau in atomic and molecular lines 
obtained with the Heterodyne Instrument for the Far Infrared on board Herschel. While CO, HjO, [Cn], and SO 
lines trace the envelope and the outflowing gas up to velocities of 33 km s~ with respect to systemic, the CN 5-4 
hyperfine structure lines at 566.7, 566.9 GHz show a narrow double-peaked profile centered at systemic velocity, 
consistent with an origin in the outer region of the compact disk of T Tau N. Disk modeling of the T Tau N disk with 
the thermo-chemical code ProDiMo produces CN line fluxes and profiles consistent with the observed ones and 
constrain the size of the gaseous disk (R oat = 110 41 -, 0 AU) and its inclination ( i = 25° ± 5°). The model indicates 
that the CN lines originate in a disk upper layer at 40-1 10 AU from the star, which is irradiated by the stellar UV 
field and heated up to temperatures of 50-700 K. With respect to previously observed CN 2-1 millimeter lines, 
the CN 5-4 lines appear to be less affected by envelope emission, due to their larger critical density and excitation 
temperature. Hence, high-/ CN lines are a unique confusion-free tracer of embedded disks, such as the disk 
of T Tau N. 
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1. INTRODUCTION 

The study of the physical and chemical structure of protoplan- 

etary disks is crucial to comprehend the formation of planetary 

systems. According to models, disks have a stratified structure, 

and hence different molecular species probe the physical and 

chemical conditions of different layers from the warm irradiated 

surface down to the cold midplane. Disks around non-embedded 
T Tauri and Herbig stars have been imaged in CO (e.g., Dutrey 

et al. 1998; Qi et al. 2004; Pietu et al. 2007), and detected in 

a number of other millimeter (mm) and submillimeter molec- 
ular lines (e.g., Dutrey et al. 1997; Thi et al. 2004). On the 

contrary, the study of disks around more embedded systems is 

obstacled by the associated envelopes and outflows which also 

emit strongly in the same molecular lines, hiding the fainter 

disk emission. 

Recent studies show that CN is a good disk tracer with 88% 

of the T Tauri disks detected in CO 2-1 which are also observed 
in the CN 2-1 lines (Oberg et al. 2010, 2011; Chapillon et al. 
2012). Moreover, Guilloteau et al. (2013) performed an IRAM 
30-m survey of T Tauri and Herbig Ae systems located mainly 
in the Taurus-Auriga region and show that the CN 2-1 lines are 
less affected by the emission from the surrounding molecular 
cloud than the lines from CO isotopologues. However, CN 2-1 


lines are still dominated by envelope/outflow emission in the 
case of actively accreting/ejecting sources, such as T Tau. 

T Tau is a multiple system driving at least two bipolar jets 
detected in optical, near-infrared forbidden lines (Bohm & Solf 
1994; Eisloffel & Mundt 1998; Solf & Bohm 1999; Herbst 
et al. 2007). The system consists of the northern component 
T Tau N (M* = 2.1 M 0 ) and of the binary system T Tau 
Sa+Sb located Oil to the south (M* = 2.1, 0.8 M 0 , separation 
Sa-Sb = O'.' 13; Kohler et al. 2008). T Tau N is optically visible 
and is surrounded by an almost face-on, intermediate mass disk 
(i < 30°, Mdi s k ~ 0.01 M 0 ) and an optically thin envelope 
(Beckwith et al. 1990; Hogerheijde et al. 1997; Akeson et al. 
1998; Ratzka et al. 2009; Guilloteau et al. 2011). Both south- 
ern components, instead, are strongly obscured (Ay — 15, 30 
mag toward T Tau Sb and Sa, respectively), possibly due to a cir- 
cumbinary envelope and, for T Tau Sa, to its almost edge-on disk 
(i ~ 72°; Duchene et al. 2005; Ratzka et al. 2009). Sa and Sb are 
not detected at mm wavelengths suggesting that their disks are 
small and of low-mass (Mdi s k ~ 10~ 5 -10~ 3 M e , R out ~ 5 AU; 
Hogerheijde et al. 1997; Ratzka et al. 2009). While the dusty 
disk of T Tau N has been mapped in the mm, there are only 
very tentative observations of its gaseous disk, e.g., in the H 2 
near-infrared lines (Gustafsson et al. 2008). Emission in the 
[Oi]63 /zm, [C 11 ] 1 57 /zm,and [Ne 11 ] 12.8 1 /imlinesis extended 
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Table 1 

Observed Lines 


Transition 

V 

(GHz) 

F up 

(K) 

HPBW 

(") 

Tpeak 

(mK) 

ypeak 
(km s _1 ) 

AV 

(km s” 1 ) 

fT mb dV 
(K km s _1 ) 

Fobs 

(W m~ 2 ) 

F mod 

(W m -2 ) 

CN 5± 9/2—7 /2 

566.731 

82 

37 

21/25 ±6 b 

6.7/9.1 ±0.5 b 

7 

0.08 ± 0.01 

5.4 ± 0.6 10~ 19 

5.7 10~ 19 

CN 5^t 11, /2— 9. /2 

566.947 

82 

37 

25/28 ± 6 b 

6.3/8.4±0.5 b 

7 

0.12 + 0.01 

8.1 ±0.6 10~ 19 

6.7 10~ 19 

0 -H 2 O 1 10 — lot 

556.936 

61 

38 

682 ±6 

8.6 ± 0.5 

52 

5.40 ± 0.02 

3.70 ± 0.02 10- 17 

3.5 10^ 19 

P-H 2 O ln-Ooo 

1113.343 

53 

19 

873 ±8 

8.6 ± 0.3 

52 

9.85 ± 0.02 

1.344 + 0.003 10~ 16 

2.2 10~ 18 

O-H 2 O 3 12 — 221 

1153.127 

249 

18 

900 ± 75 

6.2 ± 0.3 

21 

7.2 ±0.1 

1.02 + 0.02 10“ 16 

6.9 10“ 19 

CO 10-9 

1151.985 

304 

18 

13210 + 75 

7.7 ± 0.3 

38 

82.9 ± 0.2 

1.171 ±0.002 10- 15 

1.4 10~ 17 

13 CO 10-9 

1101.350 

291 

19 

1248 ± 8 

7.7 ± 0.3 

20 

4.92 ± 0.01 

6.64 ± 0.02 10- 17 

3.1 10^ 18 

SO 12i 3 — 1 1 12 

558.087 

194 

38 

19 + 6 

7.2 ± 0.5 

9 

0.08 ± 0.01 

5.3 ± 0.7 10~ 19 

7.5 10^ 20 

[C II] 2 P 3 ,2- 2 Pl/2 

1900.537 

91 

11 

1120 + 83 

7.1 ±0.2 

30 

8.4 ±0.1 

1.95 + 0.03 10~ 16 

4.8 10^ 17 


Notes. 

a Frequencies are from the Jet Propulsion Laboratory database (JPL; Pickett et al. 1998) and from the Cologne Database of Molecular Spectroscopy 
(CDMS; Muller et al. 2001) for CN. 

b For the double-peaked CN lines the intensity and velocity of both peaks is reported. 


and/or velocity-shifted, and hence likely associated with out- 
flowing gas, while molecular lines, such as high-/ CO, OH, and 
H 2 O lines, are spectrally and spatially unresolved (Spinoglio 
et al. 2000; van Boekel et al. 2009; Podio et al. 2012). Observa- 
tions of low-/ CO lines (Edwards & Snell 1982; Schuster et al. 
1993) as well as of less abundant molecular species, such as 
13 CO, C 17 0, 0 -H 2 CO, SO, and CN 2-1 are also dominated by 
envelope/outflow emission, preventing us to trace the gas in the 
disk (Guilloteau et al. 2013). 

In this Letter, we present observations of the T Tau system in 
atomic and molecular lines with higher excitation energies than 
probed by Guilloteau et al. (2013) in the mm range obtained 
with the Heterodyne Instrument for the Far Infrared (HIFI) on 
board Herschel. While CO, 13 CO, H 2 O, [Cu], and SO lines 
are still strongly dominated by emission from the surrounding 
envelope and outflows, the CN 5-4 lines allow digging out the 
faint emission from the disk of T Tau N. 

2. OBSERVATIONS AND DATA REDUCTION 

We observed T Tau (002000 = 04 h 21 m 59 ! i4, 5j2ooo = 
+19°32'06"4) with the HIFI (de Graauw et al. 2010) on board the 
Herschel Space Observatory 16 (Pilbratt et al. 2010). The obser- 
vations target the two fundamental water lines, 0 -H 2 O lio-loi 
and P-H 2 O In — Ooo in the HIFI bands 1 and 4, the 12 CO (here- 
after CO) 10-9 line in band 5, and the [Cii] 2 / > 3/2- 2 / , i/2 line 
in band 7 (OBSID: 1342249419, 1342250207, 1342249598, 
1342249647). They were acquired with a single on-source point- 
ing and in dual beam switch mode with fast chopping 3' either 
side of the target. The achieved rms noise in bands 1, 4, 5, 7 
is of 6, 8, 75, 83 mK, respectively. The Wide Band Spectrome- 
ter (WBS) and the High Resolution Spectrometer were used in 
parallel, providing a spectral resolution of 1.10 and 0.25 MHz, 
respectively. The half power beam width (HPBW) ranges from 
~11" to ~38", depending on frequency. Hence, all the three 
sources of the T Tau system (N, Sa, and Sb) are covered by 
the observations. The selected spectral settings also cover with 
the WBS the CN 5-4 / = 9/2-7/2 and / = 11/2-9/2, SO 
12i3-lli2, 13 CO 10-9, and o-H 2 0 3 i 2 -2 2 i lines. 


16 Herschel is an ESA space observatory with science instruments provided by 
European-led Principal Investigator consortia and with important participation 
from NASA. 


The HIFI data were reduced using HIPE 10. 17 Fits files 
from level 2 were then created and transformed into GILDAS 18 
format for data analysis. The spectra were baseline subtracted 
and averaged over the horizontal and vertical polarization to 
increase the signal-to-noise ratio. Antenna temperatures, 7/, 
were converted to mean beam temperature, 7 m b, using mean 
beam efficiency by Roelfsema et al. (2012). 

The properties of the detected lines (transition, frequency, 
Vo, upper level energy, E up ), the HPBW, and the measured 
parameters (peak temperature and velocity, T p ea t and V pea k, 
the full line width at la rms noise, AV, integrated intensity, 
/ T m \,dV, and line fluxes, F 0 b s 19 ) are summarized in Table 1. 

3. RESULTS FROM OBSERVATIONS 

The detected lines in Figure 1 show a variety of profiles 
indicating that they probe different gas components in the 
complex T Tau system. 

The CO and 13 CO 10-9 lines show a single-peaked profile, 
likely dominated by cloud emission, and red-shifted and blue- 
shifted wings with velocities up to —8 and +30 km s -1 probing 
outflowing gas. Assuming an average local interstellar medium 
(ISM) carbon isotope ratio of 68 ± 15 (Milam et al. 2005), 
the 12 CO/ 13 CO line ratio indicates that the CO 10-9 emission 
is optically thick, its optical depth varying between 6.7 ±1.5 
at the peak velocity down to 2.4 ± 0.8 for gas velocities of 
~5-6 km s 1 with respect to y pea k. Since the 13 CO 10-9 line is 
optically thin, it can be used to estimate the systemic velocity. 
The line peak indicates V sys = +7.7 ± 0.3 km s _1 , consistent 
with previous estimates from 13 CO and C ls O 1-0, 2-1 lines 
(Edwards & Snell 1982; Schuster et al. 1993; Guilloteau et al. 
2013). 

The [C 11 ] 2 P?,/ 2 - 2 Pi /2 and SO 12i3— 1 1 12 lines peak at 
slightly blue-shifted velocity with respect to systemic (Vp ea k = 
+7.1 ± 0.2, +7.2 ± 0.5 km s _1 ). As previously suggested by 
Guilloteau et al. (2013) SO lines trace the envelope with no 
evidence of high-velocity wings. The [C n] line instead shows 
high-velocity blue-shifted and red-shifted wings associated with 
outflowing gas. 


17 HIPE is a joint development by the Herschel Science Ground Segment 
Consortium, consisting of ESA, the NASA Herschel Science Center, and the 
HIFI, PACS, and SPIRE consortia. 

18 http://www.iram.fr/IRAMFR/GILDAS 

19 F obs = (2 K b v 3 /c 3 ) x / T mb dV x 7r(HPBW/2v4n2) 2 . 
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Figure 1. Herschel / HIFI spectra of CN 5-4 9/2 — V/2 (panel (a)), CN 5—4 
1 1/2-9/2 (panel (b)), o-H 2 0 3 i 2 -2 21 , o-H,0 1 10 -l 0 i, p-H 2 0 1 n-0 00 (panel (c)), 
CO 10-9, 13 CO 10-9 (panel (d)), SO 12 i 3 -11 12 , [Cn] 2 P 3 / 2 - 2 G /2 (panel (e)). 
Horizontal and vertical dotted lines indicate the baseline and systemic velocity 
(Vsys = +7.7 km s _1 ). The line frequency, HPBW, and upper level energy are 
labeled. For CN 5^4 the rms noise is indicated by the dashed horizontal line 
and the positions of the hyperfine structure lines by the red vertical lines. The 
line profiles predicted by the ProDiMo disk model of T Tau N are overplotted in 
black, magenta, and blue (corresponding observed profiles are in gray, magenta, 
and blue). 

(A color version of this figure is available in the online journal.) 


Podio ET AL. 

We also detect three water lines: the o-H 2 0 Iio-Iol 
0 -H 2 O 3 12—221 , and P-H 2 O I 11 -O 00 lines. The two fundamental 
water lines (£ up ~ 61, 53 K) show an absorption feature at the 
systemic velocity, due to the cloud, and a second absorption at 
~ +5.6 km s' 1 , whose origin is unclear. As already observed in 
previous works (e.g., Kristensen et al. 2012), water lines appear 
to be very sensitive to high-velocity emission showing wings 
extending up to velocities of —12 km s' 1 and +40 km s 1 . 

In contrast with previous observations of CO 6-5, 3-2, 2-1, 
and 1-0 lines (Edwards & Snell 1982; Schuster et al. 1993), the 
observed CO 10-9, 1+O, and [Cu] lines show a red- shifted wing 
which is brighter and extends to higher velocities than the blue 
one. The lines observed with HIFI have higher upper level ener- 
gies and /or critical densities than previously observed low -J CO 
lines (£ up ± 116 K), thus suggesting that the red-shifted out- 
flowing gas is denser and more excited than the blue-shifted gas. 

Our observations also cover the CN N — 5-4 hyperfine struc- 
ture lines 20 (Muller et al. 2001). We detect the three brightest / = 
912-1/2 and J = 1 1 /2-9/2 lines at 566.7 GHz and 566.9 GHz, 
while the fainter J = 9/2-9/2 components fall outside the ob- 
served range. The separation in velocity between the three de- 
tected lines is between 0.05 and 0.45 km s' 1 , and hence the 
lines are not resolved at the available resolution (0.5 km s' 1 ) 
and the observed profile is sensitive only to the kinematics of 
the emitting gas. The CN 5-4 J — 9/2-7 /2 lines show a narrow 
double-peaked profile centered at the systemic velocity with a 
total line width of ^7 km s' 1 , FWHM of 4.0 ±0.5 km s' 1 , and a 
peak separation A F sep = 2.4±0.5 km s' 1 . The profile of the J = 
11/2-9/2 blended components is also narrow and symmetric 
around the systemic velocity, even though the double-peak is not 
as clear as in the J — 9/2-7 /2 lines (see Table 1). The blue-/red- 
shifted wing shown by the J = 9/2-7 /2 and the 7=11 / 2-9/2 
components, respectively, is clearly due to the noise as it is 
below the noise level and is not detected in both components. 
Hence, the CN 5-4 profiles are consistent with emission from 
the outer region of the disk of T Tau N. The double-peak de- 
tected in the CN 5-4 lines is not seen in the 2-1 lines observed 
with the IRAM 30-m by Guilloteau et al. (2013). This is due to 
the higher excitation energy and critical density of the CN 5-4 
lines (E up ~ 82 K, n CI ~ 2 10 8 cm' 3 ), which makes them less 
affected by cloud emission, hence more effective to probe the 
disk, than CN 2-1 ( E up ~ 16 K, n CI ~ 2 10 7 cm' 3 ). Emission 
from the disks of T Tau Sa and T Tau Sb is expected to be neg- 
ligible with respect to the emission from the disk of T Tau N, as 
those disks are not detected in the mm continuum and are one to 
two orders of magnitude less massive and very small (Akeson 
et al. 1998; Ratzka et al. 2009; Guilloteau et al. 201 1). 

Assuming Keplerian rotation, a stellar mass of 2.1 M G 
and an inclination of i ~ 20°-30° (Ratzka et al. 2009), the 
peak separation of the CN lines indicates an outer disk radius 
R out (CN) < 160-350 AU. This upper limit is in agreement with 
the size of the disk as estimated from continuum maps at 1.3, 
2.7 mm (R out (dust) = 67 ± 20 AU, Guilloteau et al. 2011) and 
from the H 2 ring-like structure observed by Gustafsson et al. 
(2008) (R 0 ut(H 2 ) = 85-100 AU). 

4. MODELING CN LINES FROM THE DISK OF T Tau N 

In order to test if the CN lines originate in the disk of 
T Tau N we use a parameterized disk model calculated with the 

20 Due to the coupling of the rigid-body angular momentum, N, with the 
electronic spin, 5, and with the nuclear spin, I, angular momenta, the 
CN N = 5—4 line is split into 19 hyperfine components characterized by 
the corresponding quantum numbers J = N + S, and F = J + I . 
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Table 2 

T Tau N Disk Model: Star and Disk Parameters 


Effective temperature 

T t ff (K) 

5250 

Stellar mass 

M* (M 0 ) 

2.1 

Stellar luminosity 

U ( Lq ) 

7.3 

UV excess 

/uv 

0.1 

UV power law index 

PUV 

0.2 

X-rays luminosity 

L x (erg s -1 ) 

2 10 31 

Disk inclination 

i(°) 

25 

Disk inner radius 

Rin (AU) 

0.1 

Disk outer radius 

*out (AU) 

110 

Disk dust mass 

Miust (A^o) 

1.3 10“ 4 

Dust-to-gas ratio 

Dust-to-gas 

0.01 

Dust material mass density 

Pdust (g cm -3 ) 

2.5 

Minimum grain size 

a m m W 

0.005 

Maximum grain size 

rlmax W 

1000 

Dust size distribution index 

<7 

3.5 

Surface density X ~ r _e 

€ 

1 

Scale height at R m 

Ho (AU) 

0.0032 

Flaring index H (r) = Hq 


1.25 

Settling H(r , a) = H(r)^~ Ssa/1 

^set 

0.5 

Minimum grain size for settling 

flset (/rm) 

0.25 

Fraction of PAHs w.r.t. ISM 

/PAH 

0.01 


thermo-chemical disk modeling code ProDiMo (Woitke et al. 
2009; Kamp et al. 2010; Aresu et al. 2011). 

We adopt stellar and disk parameters as inferred from pre- 
vious studies, which well reproduce the source spectral energy 
distribution (SED; e.g., Ratzka et al. 2009). We use stellar spec- 
tral type K0 (Teff ~ 5250 K), stellar luminosity ~7.3 L 0 and 
stellar mass M* ~ 2.1 M 0 as determined by White & Ghez 
(2001). The UV spectrum (Calvet et al. 2004) is reproduced by 
an UV excess /uv = G (9 1 0—2500 A )/L* = 0.1 and a power 
law slope L , ~ a 0,2 . The effect of X-ray radiation from the 
stellar corona of T Tau N (Lx = 2 10 31 erg s -1 ; Giidel et al. 
2007) is taken into account following Aresu et al. (2011) and 
Meijerink et al. (2012). Following Ratzka et al. (2009) we as- 
sume a disk inner radius R m = 0.1 AU, a dust mixture of 
astronomical silicates (Draine & Lee 1984) and amorphous car- 
bon (Zubko et al. 1996) with relative abundances of 62.5% and 
37.5%, and a grain size distribution n(a) ~ a~ q with q = 3.5, 
where n(a) is the number of dust particles with radius a, and 
the minimum/maximum grain size are a rnm = 0.005 //in and 
a max = 1 mm. This implies a dust mass of 1.3 10~ 4 M 0 to 
reproduce the observed mm emission (Hogerheijde et al. 1997; 
Akeson et al. 1998; Guilloteau et al. 2011). Using the aver- 
age ISM dust-to-gas ratio of 0.01, the gas mass is 0.013 M 0 . 
The dust-to-gas ratio can be different from the ISM value in 
evolved disks (e.g., Thi et al. 2010; Gorti et al. 2011). However, 
as the CN lines are optically thick in the disk, their flux is in- 
dependent on the gas mass, hence on the assumed dust-to-gas 
ratio. The SED is well reproduced by assuming a disk surface 
density £ ^ r~ l and scale height H = 0.0032 AU (r/Ri n ) L25 
(Ratzka et al. 2009). To obtain a dust distribution similar to the 
two-layers model adopted by Ratzka et al. (2009), we assume 
that grains larger than a set = 0.25 //m have settled to a smaller 
scale height than the gas, H(r,d) = H (r) (a / a set )~ Sac,/1 , with 
.s'set = 0.5. As polycyclic aromatic hydrocarbon (PAH) emis- 
sion is generally not detected in T Tauri stars (e.g., Furlan et al. 
2006), the PAH fraction, /pah, is set to 0.01 with respect to the 
ISM abundance of 1 0 6 22 PAH particles/H-nucleus . Lower 
values, /pah = 10~ 3 -10~ 4 , do not affect the CN emission. The 
adopted stellar and disk parameters are summarized in Table 2. 
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2 4 6 8 10 

Line Flux (10' 19 W m‘ 2 ) 


Figure 2. Observed line properties (in red) are compared with model predictions 
(in black) for different values of the disk outer radius, S 0 ut (AU), and inclination, 
i (°), for the CN 5-4 J = 9/2-7/2 and J = 1 1 /2 9/2 lines (top and bottom 
panel, respectively). The peak separation vs. the line flux is indicated by filled 
points and solid lines connecting the model predictions for the same R out . The 
FWHM vs. the line flux is overplotted with empty circles and dotted lines. 

(A color version of this figure is available in the online journal.) 

As the CN line fluxes and profiles are very sensitive to 
the assumed disk outer radius, R out , and inclination, i, we 
proceed through two steps. First, we run a grid of models 
at low resolution (50 x 50 grid points) adopting R out and i 
values which cover the range of estimates from previous studies 
(R out = 67, 90, 100, 110, 120 AU, i = 20, 25, 30, 35, 40, 45 
AU; Stapelfeldt et al. 1998; Gustafsson et al. 2008; Ratzka et al. 
2009; Guilloteau et al. 2011; Harris et al. 2012). Then, we run 
a high resolution model (100 x 100 grid points) for R out and i 
which best fit the observations, to better resolve the chemical and 
temperature gradients and the vertical extent of the line forming 
region in the disk. The predicted FWHM and peak separation 
do not depend on the model resolution while the line fluxes are 
lower by up to ~30% when increasing the resolution. 

The line profiles and fluxes are obtained by first solving the 
statistical equilibrium with two-dimensional escape probability 
to obtain the level populations, and then using two-dimensional 
ray-tracing. As the three brightest lines of each CN 5-4 J-J' 
transition are blended and the three undetected components are 
more than one order of magnitude fainter, the model computes 
the sum of the CN 5-4 J = 9/2-7/2 and J = 11/2-9/2 lines. 
The results obtained for the grid of models (Figure 2) show 
that the CN line fluxes depend mainly on the disk size and 
increase for increasing outer radii. This is due to the fact that, 
as suggested by their profiles, the CN lines originate from the 
outer region of the disk. The best fit of the CN line fluxes 
is obtained for R out = 100-110 AU. However, as model- 
predicted fluxes are affected by ~30% uncertainty, values of 
R out between 90 and 120 AU produce CN line fluxes which 
are still in agreement with the observations. The ratio between 
the 566.9 and the 566.7 GHz lines predicted by the models, 
A mod = 1.2, is slightly lower but still in agreement with the 
observed one (R 0 bs = 1.4 ± 0.2). The peak separation and 
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1 10 100 
r [AU] 


Figure 3. Disk region from which 50% of the CN 5 — 4- y = 9/2 — 7 /2 at 566.7 GHz 
(in orange), and the CN 5-4 7=11 /2-9/2 at 566.9 GHz (in red) line emission 
arises in the disk model of T Tau N. Note that the emitting regions of the two CN 
lines overlap almost completely. The gray color indicates the CN density, «cn 
( cm -3 ), the dotted black and blue curves the gas temperature and total hydrogen 
number density «(h) = «h + 2;!h 2 . 

(A color version of this figure is available in the online journal.) 

FWHM predicted by the models depend on both R tml and i. 
For R out — 90-120 AU the observed AF sep and FWHM are 
reproduced if the disk inclination is low (i = 20° — 30°). This is 
in agreement with the conclusions by Akeson et al. (2002) and 
Ratzka et al. (2009) based on the modeling of the SED and near-, 
mid-infrared visibilities. Larger inclination angles of 40° -45° 
(Stapelfeldt et al. 1998; Guilloteau et al. 2011) are excluded 
based on the observed CN profiles. 

Following the results obtained from the grid of models, we 
adopt R lml =110 AU, i = 25° and compute the fluxes and 
profiles of all the lines covered by our observations using the 
high resolution model. Figure 3 shows the region in the disk 
where 50% of the CN N = 5-4 / = 9/2-7/2 and / = 11/ 
2-9/2 lines originate according to our model. This is obtained 
using vertical escape probability and without accounting for disk 
inclination. The model indicates that the CN lines are excited 
in a disk upper layer located at 40-1 10 AU distance from the 
star, which is irradiated by the stellar UV field. This heats the 
gas up to temperatures of ~ 5 0-7 00 K, while the gas density is 
~3 10 7 — 3 I 0 9 cm” 3 , thus the CN lines are almost thermalized 
and optically thick (r ~ 10 2 — 10 3 ) in the line emitting region. 
In this disk layer, CN is produced by H + CN + , and C + 
NO reactions, and photo-dissociation of HCN, while the main 
destruction route is photo-dissociation of CN. 

In Figure 1 and Table 1, the model-predicted line profiles and 
fluxes are compared with the observed ones. For CN lines, the 
model predicts an FWHM of 3.9 km s” 1 and peak separation 
of 2.6 km s” 1 , in agreement with observations (FWHM = 4.0, 
4.2 ± 0.5 km s” 1 , Ay sep = 2.4, 2.1 ± 0.5 km s” 1 ). Also the line 
fluxes agree with observations within a factor ~0.9, 1.2. The 
sum of the fluxes of the N — 2-1 components, instead, is about 
an order of magnitude lower than what observed by Guilloteau 
et al. (2013), who suggest that the line is dominated by envelope 
emission. The CO, 13 CO, fLO, and SO line fluxes predicted by 
the disk model are from a factor four to two orders of magnitude 
smaller than observed, as well as for the high-/ CO, H 2 O, OH, 
and atomic [Oi], [C 11 ] lines previously observed with PACS 
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(see Table 4 by Podio et al. 2012). This is a further evidence that 
these lines are dominated by envelope and outflow emission, 
as already suggested by the fact that [Oi], [C 11 ] emission is 
spatially extended with PACS Q^9") and by the broad profiles 
of CO, H 2 O, SO, and [Cii] obtained with HIFI. 

5. CONCLUSIONS 

The Herscliel / HIFI observations of the T Tau system show 
emission in a number of molecular and atomic lines. The origin 
of the emission lines in embedded, accreting/ejecting sources 
is highly debated, being crucial to identify a tracer to dig out the 
faint disk emission (e.g., Podio et al. 2012, 2013). In the case of 
T Tau, the CO, H 2 O, and [Cii] lines clearly trace high-velocity 
outflowing gas. By contrast with those lines and with previously 
observed CN 2-1 lines (Guilloteau et al. 2013), the CN 5-4 lines 
show narrow double-peaked profiles centered at the systemic 
velocity, suggesting an origin in the outer disk of T Tau N. Disk 
modeling predicts CN line fluxes and profiles in agreement with 
observed ones and constrains the size of the gaseous disk of 
T Tau N (R out = 1 1 0//, 1 , AU) and its inclination (i = 25° ± 5°). 
This study demonstrates that high-/ CN lines are a unique 
tool to probe the gaseous disk of strongly accreting/ejecting 
sources and paves the way for future observations of embedded 
disks with ALMA. 

L.P. and S.D.B. acknowledge funding from the European FP7 
(PIEF-GA-2009-253896) and the National Science Foundation 
(AST-0954811). 
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